SHOCK WAVES AND COMPACTONS FOR 
FIFTH-ORDER NONLINEAR DISPERSION EQUATIONS 

VICTOR A. GALAKTIONOV 

Abstract. The following question is posed: to justify that the standing shock wave 

„ / N . f — 1 for x < 0, 

■J-(x) = — signx = — < . n 
y ' b | 1 for i > 0, 

is a correct "entropy" solution of fifth-order nonlinear dispersion equations (NDEs), 

u t = — (uu x ) xxxx and u t = - {uu xxxx ) x in R x R+. 

These two quasilinear degenerate PDEs are chosen as typical representatives, so other 
similar (2m + l)th-order NDEs with no divergence structure admit such shocks. 

As a related second problem, the opposite shock S+{x) = —S-(x) = sign a; is shown 
to be a non-entropy solution that gives rise to a continuous rarefaction wave for t > 0. 

Formation of shocks is also studied for the fifth-order in time NDE 

uttttt = (uu x ) xxxx • 

On the other hand, related NDEs are shown to admit smooth compactons, e.g., for 
u t = -{\u\u x ) xxxx + \u\u x in RxM + , 
which are of changing sign. Nonnegative ones are nonexistent in general (not robust). 



1. Introduction: nonlinear dispersion PDEs, Riemann's Problems, and 

main directions of study 

1.1. Three main problems: entropy shocks and rarefaction waves for fifth-order 
NDEs. Let us introduce our basic models that are five fifth-order nonlinear dispersion 
equations (NDEs). These are ordered by numbers of derivatives inside and outside the 
quadratic differential operators involved: 



(1.1) u t = -uu xxxxx (NDE-(5,0)), 

(1.2) u t = -{uu xxxx ) x (NDE-(4,1)), 

(1.3) u t = -{uu xxx ) xx (NDE-(3,2)), 

(1.4) u t = -{uu xx ) xxx (NDE-(2,3)), 

(1.5) u t = -{uu x ) xxxx (NDE-(1,4)). 
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The only fully divergent operator is in the last NDE-(1,4) that, being written as 

(1.6) ut = ~{uu x ) xxxx = -\ (u 2 ) xxxxx (NDE-(1,4) = NDE-(0,5)), 

becomes also the NDE-(0,5), or simply the NDE-5. This completes the list of such 
quasilinear degenerate PDEs under consideration. 

Before explaining the physical significance of such models and their role in general PDE 
theory, we pose those three main problems for the above NDEs: 

(I) Problem "Blow-up to S-" (Section |2j): to show that the shock of the shape —sign a; 
can be obtained by blow-up limit from regular (at least, continuous) solutions U-(x,t) of 
(TP]) - (TO]) mix (0,T), so that 

(1.7) u-(x,t) ^£_(x) = -signx= ^ x>^0, aS t ~^ T ' m L iocW- 

(II) Riemann's Problem S + (RP+) (Sectional): to show that the shock 



^1.8) S+(x) = sign x 



-1 for x < 0, 
1 for x > 0, 



for NDEs (li.Jjl generates "rarefaction waves" that are continuous for t > 0. 

(III) Riemann's Problem SL (RP-) (Section HJ): to show that 

(1.9) S-{ x ) is an "entropy" shock wave, and S+(x) is not. 

(IV) Problem: nonuniqueness/ entropy (Section Ej): to show that a single point "gradient 
catastrophe" for the NDE ( j 1 . 51) leads to a shock wave, which is principally nonunique. 

In Section [6], we discuss first these problems in application to various NDEs including 
the following rather unusual one: 

(1-10) Uuttt — (uu x ) xxxx , 

which indeed can be reduced to a first-order system that, nevertheless, is not hyperbolic. 
The main convenient mathematical feature of (11.101) is that it is in the normal form, so 
it obeys the Cauchy-Kovalevskaya theorem that guarantees local existence of a unique 
analytic solution (this adds extra flavour to our 5-entropy test). Regardless this, (jl.lOp is 
shown to create shocks in finite time and rarefaction waves for other discontinuous data. 

In Section [7J as the last "opposite to shocks problem" and as a typical example, we 
consider the following perturbed version of the NDE (11. 5ft : 

(1-11) u t = -{\u\u x ) xxxx + \u\u x , 

which is written for solutions of changing sign by replacing u 2 by the monotone function 
\u\u. This is essential, and we show that (11.111) admits compactly supported travelling 
wave (TW) solutions of changing sign called compactons. More standard in literature 
nonnegative compactons of fifth-order NDEs such as (11.111) and others are shown to be 
nonexistent in general in the sense that these are not robust, i.e., do not exhibit continuous 
dependence on parameters of PDEs and/or small perturbations of nonlinearities. 



1.2. A link to classic entropy shocks for conservation laws. The above problems 
(I) (III) are classic for entropy theory of ID conservation laws from the 1950s. Shock 
waves first appeared in gas dynamics that led to mathematical theory of entropy solutions 
of the first-order conservation laws such as Euler's equation 

(1.12) u t + uu x = in R x R + . 

Entropy theory for PDEs such as (11.121) was created by Oleinik [3"0"l ETJ and Kruzhkov [2B] 
(equations in M. N ) in the 1950-60s; see details on the history, main results, and modern 
developments in the well-known monographs [H [9j 02]. Note that first analysis of the 
formation of shocks was performed by Riemann in 1858 [51] ; see the history in [3j. 

According to entropy theory for conservation laws such as (I1.12p . it is well-known that 
(11. 9p holds. This means that 

(1.13) U-(x, t) = S-(x) = — sign x 

is the unique entropy solution of the PDE (I1.12p with the same initial data S~(x). On 
the contrary, taking S^-type initial data yields the continuous rarefaction wave with a 
simple similarity piece- wise linear structure, 

{—1 for x < —t, 
j for \x\ < t, 
1 for x > t. 

Our goal is to justify the same conclusions for fifth-order NDEs, where of course the 
rarefaction wave in the RP+ is supposed to be different from that in (I1.14p . 

We now return to main applications of the NDEs. 

1.3. NDEs from theory of integrable PDEs and water waves. Talking about odd- 
order PDEs under consideration one should mention that these naturally appear in classic 
theory of integrable PDEs, with those representatives as the KdV equation, 

(1.15) u t + uu x = u xxx , 
the fifth- order KdV equation, 

u t + u xxxxx + 30 u 2 u x + 20 u x u xx + 10 uu xxx = 0, 

and others from shallow water theory. These are semilinear dispersion equations, which 
being endowed with smooth semigroups, generate smooth flows, so discontinuous weak 
solutions are unlikely; see references in [201 Ch. 4]. 

The situation is changed for the quasilinear case. In particular, for the quasilinear 
Harry Dym equation 

(1.16) u t = u 3 u xxx , 

which is one of the most exotic integrable soliton equations; see [20| § 4.7] for survey and 
references therein. Here, (11.161) indeed belongs to the NDE family, though it seems that 
semigroups of its discontinuous solutions were never under scrutiny. In addition, integrable 
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equation theory produced various hierarchies of quasilinear higher-order NDEs, such as 
the fifth-order Kawamoto equation [25] 

(1.17) Ut U Uxxxxx ~\~ 5 U U x U X xxx 10 It U XX U XXX . 

We can extend this list talking about possible quasilinear extensions of the integrable 

Lax's seventh-order KdV equation 

u t + [35u 4 + 70(u 2 u xx + u{u x ) 2 ) + 7(2uu xxxx + 3(u xx ) 2 + Au x u xxx ) + u xxxxxx ] x = 0, 

and the seventh-order Sawada-Kotara equation 

u t + [63u 4 + 63(2u 2 u xx + u{u x ) 2 ) + 1\[uu xxxx + [u xx ) 2 + u x u xxx ) + u xxxxxx ] x = 0; 

see references in [201 P- 234]. 

Modern mathematical theory of odd-order quasilinear PDEs is partially originated and 
continues to be strongly connected with the class of integrable equations. Special ad- 
vantages of integrability by using the inverse scattering transform method, Lax pairs, 
Liouville transformations, and other explicit algebraic manipulations made it possible to 
create rather complete theory for some of these difficult quasilinear PDEs. Nowadays, 
well-developed theory and most of rigorous results on existence, uniqueness, and various 
singularity and non-differentiability properties are associated with NDE-type integrable 
models such as Fuchssteiner-Fokas-Camassa-Holm (FFCH) equation 

(1.18) (I - D 2 )u t = -3uu x + 2u x u xx + uu xxx = -(I - D 2 x )(uu x ) - [u 2 + \ [u x ) 2 \ . 

Equation (11.181) is an asymptotic model describing the wave dynamics at the free surface 
of fluids under gravity. It is derived from Euler equations for inviscid fluids under the 
long wave asymptotics of shallow water behavior (where the function u is the height of 
the water above a flat bottom). Applying to fl!.18[) the integral operator (/ — D 2 )~ l with 
the L 2 -kernel u(s) = ^e~' s ' > 0, reduces it, for a class of solutions, to the conservation 
law (11.121) with a compact first-order perturbation, 

(1.19) u t + uu x = -[u* (u 2 + \{u x ) 2 )] x . 

Almost all mathematical results (including entropy inequalities and Oleinik's condition 
(E)) have been obtained by using this integral representation of the FFCH equation; see 
a long list of references given in [20, p. 232]. 

There is another integrable PDE from the family with third-order quadratic operators, 

(1.20) u t - u xxt = auu x + (3u x u xx + uu xxx (a, (3 G R), 

where a = — 3 and (3 = 2 yields the FFCH equation fll . 18j) . This is the Degasperis-Procesi 
equation for a = —4 and (3 = 3, 

(1.21) u t - u xxt = -Auu x + 3u x u xx + uu xxx . 

On existence, uniqueness (of entropy solutions in L 1 n BV), parabolic e-regularization, 
Oleinik's entropy estimate, and generalized PDEs, see [5]. Besides (11.181) and (11.211) . the 
family fl 1 . 2 [) does not contain other integrable entries. A list of more applied papers 
related to various NDEs is also available in [201 Ch. 4]. 
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1.4. NDEs from compacton theory. Other important application of odd-order PDEs 
are associated with compacton phenomena for more general non-integrable models. For 
instance, the Rosenau-Hyman (RH) equation 

(1.22) ut = {u% xx + {u 2 ) x 

which has special important application as a widely used model of the effects of nonlinear 
dispersion in the pattern formation in liquid drops [39]. It is the K(2,2) equation from 
the general K(m, n) family of the following NDEs: 

(1.23) ut = {u n ) xxx + {u m ) x (u>0), 

that describe phenomena of compact pattern formation, [351 EH] . Such PDEs also appear 
in curve motion and shortening flows [38J. Similar to well-known parabolic models of 
porous medium type, the K(m, n) equation (11 .23[) with n > 1 is degenerated at u — 0, 
and therefore may exhibit finite speed of propagation and admit solutions with finite 
interfaces. The crucial advantage of the RH equation (11. 22ft is that it possesses explicit 
moving compactly supported soliton-type solutions, called compactons [39], which are 
travelling wave (TW) solutions to be discussed for the PDEs under consideration. 

Various families of quasilinear third-order KdV-type equations can be found in [4], 
where further references concerning such PDEs and their exact solutions can be found. 
Higher-order generalized KdV equations are of increasing interest; see e.g., the quintic 
KdV equation in [23] and [17], where the seventh-order PDEs are studied. For the K(2, 2) 
equation (11.221) . the compacton solutions were constructed in [35] . 

More general B(m, k) equations, coinciding with the K(m, k) after scaling, u t +a(u m ) x = 
fi>(u k ) xxx also admit simple semi-compacton solutions [40J, as well as the Kq(m,u) non- 
linear dispersion equation (another nonlinear extension of the KdV) [SS] 

u t + (u m ) x + [u 1 - u {u u u x ) x ] x = 0. 

Setting m = 2 and lu — | yields a typical quadratic PDE u t + {u 2 ) x + uu xxx + 2u x u xx = 
possessing solutions on standard invariant trigonometric-exponential subspaces, where 
u(x, t) = Co(t) + Ci(t)cosXx + C2(t)sm\x and {Cq,C\,C2} solve a nonlinear 3D dy- 
namical system. Combining the K(m,n) and B(m,k) equations gives the dispersive- 
dissipativity entity DD(k, m, n) [37] u t + a{u m ) x + (u n ) xxx = n(u k ) xx that can also admit 
solutions on invariant subspaces for some values of parameters. 
For the fifth-order NDEs, such as 

(1.24) u t = a(u 2 ) xxxxx + f3{u 2 ) xxx + ^{u 2 ) x in K x E + , 

compacton solutions were first constructed in [10] . where the more general K(m,n,p) 
family of PDEs, Ut+f3i(u m ) x +f32(u n ) xxx +f33D^(u p ) = 0, with m, n, p > 1, was introduced. 
Some of these equations will be treated later on. Equation (11.241) is also associated with 
the family Q(l,m,n) of more general quintic evolution PDEs with nonlinear dispersion, 

(1.25) ut + a(u m+1 ) x + u [u{u n ) xx ] % + 5 [u{u l ) xxxx ] x = 0, 

possessing multi-hump, compact solitary solutions [41] . 
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Concerning higher-order in time quasilinear PDEs, let us mention a generalization of 
the combined dissipative double- dispersive (CDDD) equation (see, e.g., [33J) 

(1-26) u tt = au xxxx + f3u xxtt + j(u 2 ) xxxxt + 5{u 2 ) xxt + e(u 2 ) t , 

and also the nonlinear modified dispersive Klein-Gordon equation (mKG(l,n, k)), 

(1.27) u tt + a(u n ) xx + b(u k ) xxxx = 0, n, k > 1 (u > 0); 

see some exact TW solutions in [24J. For b > 0, (11.271) is of hyperbolic (or Boussinesq) type 
in the class of nonnegative solutions. We also mention related 2D dispersive Boussinesq 
equations denoted by B(m,n, k,p) [46J, 

(u m ) tt + a(u n ) xx + (3{u k ) xxxx + l(u p ) yyyy = in K 2 xi 

See [20j Ch. 4-6] for more references and examples of exact solutions on invariant subspaces 
of NDEs of various types and orders. 

1.5. On canonical third-order NDEs. Until recently, quite a little was known about 
proper mathematics concerning discontinuous solutions, rarefaction waves, and entropy 
approaches, even for the simplest third-order NDEs such as (11.221) or (see |19j ) 

(1.28) u t = (uu x ) xx . 

However, the smoothing results for sufficiently regular solutions of linear and nonlinear 
third-order PDEs are well know from the 1980-90s. For instance, infinite smoothing results 
were proved in [6] (see also |22j ) for the general linear equation 

(1.29) u t + a(x,t)u xxx = (a(x,t) > c> 0), 
and in [7] for the corresponding fully nonlinear PDE 

(1.30) u t + f{u 

XXX) 'U'XXI ^X1 ^1 •£) 

t) = {f Uxxx >c>0); 

see also [2] for semilinear equations. Namely, for a class of such equation, it is shown that, 
for data with minimal regularity and sufficient (say, exponential) decay at infinity, there 
exists a unique solution u(x,t) 6 C£° for small t > 0. Similar smoothing local in time 
results for unique solutions are available for equations in M 2 , 

(1.31) u t + f(D 3 u, D 2 u, Du, u, x, y, t) = 0; 

see [27] and further references therein. 

These smoothing results have been used in [16] for developing some 5-entropy concepts 
for discontinuous solutions by using techniques of smooth deformations. We will follow 
these ideas applied now to shock and compacton solutions of higher-order NDEs and 
others. 

2. (I) Problem "Blow-up": EXISTENCE OF SIMILARITY SOLUTIONS 

We now show that Problem (I) on blowing up to the shock S-(x) can be solved in 
a unified manner by constructing self-similar solutions. As often happens in nonlinear 
evolution PDEs, the refined structure of such bounded and generic shocks is described in 
a scaling-invariant manner. 
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2.1. Finite time blow-up formation of the shock wave S~(x). One can see that all 
five NDEs fll.ip - fll.5l) admit the following similarity substitution: 

(2.1) u-{x,t) = g(z), z = x/(-t)*, 

where, by translation, the blow-up time in reduces to T = 0. Substituting (12. II) into the 
NDEs yields for g the following ODEs in M, respectively: 



(2.2) 


gg® = -\sfz, 


(2.3) 


{gg^)' = -\g'z 


(2.4) 


(g 9 "')" = -lg'z, 


(2.5) 


{gg'T = -\9'z, 


(2.6) 


{ggT = -\9'z 



with the following conditions at infinity for the shocks S-\ 
(2.7) 0(=foo) = ±1. 

In view of the symmetry of the ODEs, 



(2i 



9 >-> -g, 

Z I — > —z. 



it suffices to get odd solutions for z < posing anti-symmetry conditions at the origin, 
(2.9) £(0)=/(0)=s (4) (0) = 0. 

2.2. Shock similarity profiles exist and are unique: numerical results. Before 
performing a rigorous approach to Problem (I), it is convenient and inspiring to check 
whether the shock similarity profiles g(z) announced in (12. ip actually exist and are unique 
for each of the ODEs fl2.2l) - fl2.6p . This is done by numerical methods that supply us 
with positive and convincing conclusions. Moreover, these numerics clarify some crucial 
properties of profiles that will clarify the actual strategy of rigorous study. 

A typical structure of this shock similarity profile g(z) satisfying (12. 2p . (12. 9p is shown 
in Figure [TJ As a key feature, we observe a highly oscillatory behaviour of g(z) about ±1 
as z — > =Foo, that can essentially affect the topology of the announced convergence (11 .7j) . 
Therefore, we will need to describe this oscillatory behaviour in detail. In Figure [2], we 
show the same profile g(z) for smaller z. It is crucial that, in all numerical experiments, 
we obtained the same profile that indicates that it is the unique solution of (12. 2p . (12. 9p . 

Figure [3](a)-(d) show the shock similarity profiles for the rest of NDEs fll.2p - fll.5p . They 
differ from each other rather slightly. 

Remark: on regularization in numerical methods. For the fifth-order NDEs, this 
and further numerical constructions are performed by MatLab by using the bvp4c solver. 
Typically, we take the relative and absolute tolerances 



(2.10) 



Tols = 10~ 4 . 
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Shock similarity profile for NDE-(5,0): u=-u u 

3 r I xxxxx 



gg (5) =-g'z/5 
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FIGURE 1. The shock similarity profile g(z) as the unique solution of the prob- 
lem ([22D, CUD; z e [-200,200]. 



Shock similarity profile for NDE-(5,0): u=- u u 

t xxxxx 




FIGURE 2. The shock similarity profile g(z) as the unique solution of the prob- 
lem flZ2]), (JZS]); z £ [-20,20]. 



Instead of the degenerate ODE (12.21) (or others), we solve the regularized equation 



(2.11) 
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(5) 



. (\g'z). 



with the regularization parameter v — 10 



where the choice of small v is coherent with the tolerances in (12.101) . Sometimes, we will 
need to use the enhanced parameters Tols ~ v ~ 10~ 7 or even 10 -9 . 
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Shock similarity profile for NDE-(4,1 ): u t =-(u Shock similarity profile for NDE-(3,2): u t =- (u u xxx ). 




(c) equation (|2.5p (d) equation (|2.6p 

FIGURE 3. Shock similarity profiles as solutions of (|2.3p - (|2.6p . (|2.9p respectively. 
For comparison, dotted lines denote the profile from Figures Q] and [2j 



2.3. Justification of oscillatory behaviour about equilibria ±1 and other asymp- 
totics. Thus, the shock profiles g(z) are oscillatory about +1 as z — > — oo. In order to 
describe these oscillations in detail, we linearize all the ODEs (I2.2l) - (l2.6p about the regular 
equilibrium to get the linear ODE 

(2.12) g® = -yz. 

This equation reminds that for the rescaled kernel F(z) of the fundamental solution of 
the corresponding linear dispersion equation, 

(2.13) u t = —u xxxxx in R x R + . 

The fundamental solution of the corresponding linear operator + D x in (I2.13P has the 
standard similarity form 

(2.14) b(x,t) =t~*F(y), with y = x/t 1 ^ 5 , 
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where F(y) is a unique solution of the ODE problem 

(2.15) B 5 F = -F (5) + | (Fy)' = in R, / F = 1, or = | on integration. 

The precise asymptotics of small solutions of ( 12.121) as z —>■ — oo is as follows: g(z) ~ 
e a l z l 5 4 , where a 4 = Choosing the purely imaginary root with Re a = gives a more 
refined WKBJ-type asymptotics, 

(2.16) g(z) = 1 + Co|z|~^ [Asin(ao|^|3) + B cos(ao|^|3)] + ... . 
This asymptotic behaviour implies two conclusions: 

(2.17) g(z)-l^L 1 (R_) 

and that the total variation of g(z) (and hence of u-(x, t) for any t < 0) is infinite. Setting 
|^| 5 / 4 = v in the integral yields 

+oo oo oo 

(2.18) I^(')Itv = / |0'( 2 )|cLz~ / z ~* I cosz^l dz ~ J df = oo. 

— oo 

This is in striking contrast with the case of conservation laws ( 11 . 121) . where finite total 
variation approaches and Helly's second theorem (compact embedding of sets of bounded 
functions of bounded total variations into L°°) used to be key; see Oleinik's pioneering 
approach [30]. In view of the presented properties of the similarity profile g{z), the 
convergence in (jl.7p takes place for any igK, uniformly in R\ (— //, //),//> small, and 
in L^ oc (R) for p G [1, oo), that, for convenience, we fix in the following: 

Proposition 2.1. For the shock similarity profile g(z) the convergence ?p with T = 0: 

(i) does not hold in L 1 (R), and 

(ii) does hold in L loc (R) ; and moreover, for any fixed finite I > 0, 

(2.19) |M-,t)-S_(-)IUi(- M) =O((-^)^0 as t->0- 

Proof of (12.191) is the same as in (12.181) with a finite interval of integration: for 1 = 1, 

II ' || z,i(— 1,1) ~ (— £) 5 / z s | cosz4 1 dz ~ (—t) s f v io | cosf | dz ~ (— t) « . 

Finally, note that each g(z) has a regular asymptotic expansion near the origin. For 
instance, for the first ODE (12.21) . there exist solutions such that 

(2.20) g(z) = Cz + Dz 3 -^- z 5 + ^z 7 + ..., 

where C < and D e R are some constants. The uniqueness of such asymptotics is traced 
out by using Banach's Contraction Principle applied to the equivalent integral equation 
in the metric of C(—fi,fi), with fi > small. Moreover, it can be shown that (12.201) is 
the expansion of an analytic function. Other ODEs admit similar local representations of 
solutions. 

In addition, we need the following scaling invariance of the ODEs (I2.2l) - ()2.6p : if g\{z) 
is a solution, then 

(2.21) 9a(z) = fl5 5 , i(f) is a solution for any a 7^ 0. 
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2.4. Existence of shock similarity profiles. Using the asymptotics derived above, we 
now in a position to prove the following: 



Proposition 2.2. The problem i \2.7[ , ( \2.9h for ODEs ( TO) -( TO) admits a solution g(z), 



which is an odd analytic function. 

Notice that Figures above clearly convince that, moreover, 

(2.22) g(z) > for z < 0, 

which is difficult to prove rigorously; see further comments below. Actually, f 1 2 . 2 2 j) is not 
that important for the key convergence (11.71) . 

Proof. As above, we consider the first ODE (12. 2p only. We use the shooting argument 
using the 2D bundle of asymptotics (12.201) . By scaling (12.211) . we put C — — 1, so, actually, 
we deal with the one-parameter shooting problem with the ID family of orbits satisfying 

(2.23) g(z; D) = -z + D z 3 - ± z 5 - JjL z 7 + ... , Del. 
It is not hard to check that, besides constant unstable equilibria, 

(2.24) g(z) -> C- > as z -> -oo, 
the ODE (12. 3p admits an unbounded stable behaviour given by 

(2.25) g(z) ~ g*(z) = — z 5 — >■ +oo as z — > — oo. 

This determines the strategy of the ID shooting via the family (12.231) : 

(i) obviously, for all D <C —1, we have that g(z;D) > is monotone decreasing and 
approaches the stable behaviour (I2.25p . and 

(ii) on the contrary, for all D ^> 1, g{z\ D) gets non-monotone and has a zero for some 
finite z = z (D) < 0, z (D) — > _ as D — >• — oo, and eventually approaches (12.251) . but 
in an essentially non-monotone way. 

It follows from different and opposite "topologies" of the behaviour announced in (i) 
and (ii) that there exists a constant D such that g(z;D Q ) does not belong to those two 
sets of orbits (both are open) and hence does not approach g*(z) as z —>■ — oo at all. This 
is precisely the necessary shock similarity profile with some unknown properties. □ 

This ID shooting approach is explained in Figure H] obtained numerically, where 

(2.26) D = 0.069192424... . 

It seems that as D — » Dq , the zero of g(z; D) must disappear at infinity, i.e., 

(2.27) z {D) -> -oo as D -> D+, 

and this actually happens as Figure H] shows. Then this would justify the positivity 
(12.221) . Unfortunately, in general (i.e., for similar ODEs with different sufficiently arbitrary 
nonlinearities), this is not true, i.e., cannot be guaranteed by a topological argument. So 
that the actual operator structure of the ODEs should be involved, so, theoretically, the 
positivity is difficult to guarantee in general. Note again that, if the shock similarity 
profile g(z) would have a few zeros for z < 0, this would not affect the crucial convergence 
property such as (11. 7\i . 
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Shooting shock similarity profile for C=-1 : D =0. 0691 92424... 

2 
1.8 
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0.4 
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-20 -15 -10 -5 

z 

FIGURE 4. Shooting the shock similarity profile g(z) via the family (f!T23j) ; D = 0.069192424.... 

2.5. Self-similar formation of other shocks. NDE-(1,4). Let us first briefly consider 
the last ODE ( 12.61) for the fully divergent NDE (11.51) . Similarly, by the same arguments, 
we show that, according to (12.11) . there exist other shocks as non-symmetric step-like 
functions, so that, as t — > 0~, 

( C_ > for x < 0, 

(2.28) u_(x,t)^<C for x = 0, 

(C+ < for x > 0, 

where C_ ^ —C + and C ^ 0. Figure [5] shows a few of such similarity profiles g(z), where 
three of these are strictly positive. The most interesting is the boldface one with C_ = 1.4 
and C + = that has the finite right-hand interface at z = z pa 5, with the expansion 

(2.29) g{z) = - _ z f + (i + (i)) ^ " as Z -> z - 

It follows that this g(z) < near the interface so the function changes sign there, which 
is also seen in Figure [5] by carefully checking the shape of profiles above the boldface one 
with the finite interface bearing in mind a natural continuous dependence on parameters. 

NDE-(5,0). Consider next the first ODE (Q for the fully non-divergent NDE-(5, 0) 
(II. ip . We again can describe formation of shocks (12.281) ; see Figure [61 The boldface 
profile with C_ = 1.4 and C + = has finite right-hand interface at z = z$ ~ 5, with a 
different expansion 

(2.30) g{z) = § f(z - z) 4 \ \n(z - z)\(l + o(l)) -> 0+ as z -> z ". 
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Shock similarity profiles for NDE-(1 ,4)=(0,5): u=- (u u x ). 
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FIGURE 5. Various shock similarity profiles g(z) as solutions of the problem 
(USD, ([IS]). 



Shock similarity profiles for NDE-(5,0): u=- u u 




z 

FIGURE 6. Various shock similarity profiles g(z) as solutions of the problem 

CL2D, CS1. 
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Shock similarity profiles for u =- ((1+u )u 




FIGURE 7. Various shock similarity profiles g(z) satisfying the ODE (|2.32j) . 

2.6. On shock formation for a uniform NDE. Here, as a key example (to be contin- 
ued), we show shocks for uniform NDEs such as the fully divergent one 

(2.31) u t = -((1 + U 2 )u x )xxxx- 

This equation is non-degenerate, so represents a "uniformly dispersive" equation. The 
ODE for self-similar solutions (12. ip then takes the form 

(2.32) (( l + ^ y) (4) = _i^_ 

The mathematics of such equations is similar to that in Section [2741 In Figured we present 
a few shock similarity profiles for (I2.32p . Note that both shocks S±(x) are admissible, 
since for the ODE (and the NDE (I2.3ip ). we have, instead of symmetry d 2 . 8 j) . 

(2.33) — g(z) is also a solution. 



3. Riemann Problem S + : SIMILARITY RAREFACTION WAVES 

Using the reflection symmetry of all the NDEs under consideration, 

fit i — ► — u, 

3.1 <^ 

we conclude that these admit global similarity solutions defined for all t > 0, 
(3.2) u+(x, t) = g(z), with z = x/t$. 

Then g(z) solves the ODEs (I2.2l) - (l2.6p with the opposite terms 
(3.3) 
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on the right-hand side. The conditions (12. 7p also take the opposite form 

(3.4) /(±oo) = ±1. 

Thus, these profiles are obtained from the blow-up ones in (12.11) by reflection, 

(3.5) if g(z) is a blow-up shock profile in (I2.ip . then g{—z) is rarefaction one in (13. 2p . 

These are sufficiently regular similarity solutions of NDEs that have the necessary initial 
data: by Proposition I2.1( ii). in L\ oc , 

(3.6) u + (x, t) — > S+(x) as t — > + . 

Other profiles g(—z) from shock wave similarity patterns generate further rarefaction 
solutions including those with finite left-hand interfaces. 

4. Riemann Problem £>_: TOWARDS 5-ENTROPY TEST 

4.1. Uniform NDEs. In this section, for definiteness, we consider the fully non-divergent 
NDE (HID, 

(4.1) ut = A(u) = -uu xxxxx in Kx (0, T), u(x, 0) = u (x) G C °°(R). 

In order to concentrate on shocks and to avoid difficulties with finite interfaces or transver- 
sal zeros at which u = (these are weak discontinuities via non- uniformity of the PDE), 
we deal with strictly positive solutions, say, 

(4.2) ^ < u < C, where C > 1 is a constant. 

Remark: uniformly non-degenerate NDEs. Alternatively, in order to avoid the 
assumptions like (14. 2p . we can consider the uniform equations such as, e.g., 

(4.3) u t = -(1 + u 2 )u xxxxx , 

for which no finite interfaces are available. Of course, (14. 3 j) admits similar blow-up simi- 
larity formation of shocks by (12.11) . In Figure [U we show a few profiles satisfying 

(4.4) (i + <rV 5) = -i^> zeR - 

Recall that, for ( 14.41) . (12.331) holds, so both S±(x) are admissible and entropy (see below). 

4.2. On uniqueness, continuous dependence, and a priori bounds for smooth 
solutions. Actually, in our 5-entropy construction, we will need just a local semigroup 
of smooth solutions that is continuous is L\ oc . The fact that such results are true for 
fifth-order (or other odd-order NDEs) is easy to illustrated as follows. 

One can see that, since (14.11) is a dispersive equation, which contains no dissipative 
terms, the uniqueness follows as for parabolic equations such as 

Ut = —uu xxxx or u t = uu xxxxxx (in the class {^<w<C}). 

Thus, we assume that u(x, t) solves (14. ip with initial data uo(x) G H W (R), satisfies 
(Q|) . and is sufficiently smooth, u G L°°([0, T], H 10 (R)), u t G L°°([0, T], if 5 (M)), etc. 
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Shook similarity profiles for u=- (1+u )u 

1 r t v ' xxxxx 




-100 -80 -60 -40 



FIGURE 8. Various shock similarity profiles g{z) satisfying the ODE (14. 4 



Assuming that v(x, t) is the second smooth solution, we subtract equations to obtain for 
the difference w = u — v the PDE 

(4-5) Wi uw xxxxx v xxxxx w. 

We next divide by u > ^ > and multiply by w in L 2 , so integrating by parts that 
vanishes the dispersive term w xxxxx yields 



(4.6) 



r ww t — 1 _d_ f I 1 f ut 
J u ~ 2 At J u ~ t ~ 2 J u 2 



I 



w 



Therefore, using f!4.2[) and the assumed regularity yields 



(4.7) 



1 d_ 

2 dt 



[ = r ( _ i _ fe^) w 2 < c J h! 

Jm JV 2m^ u / — 



where the derivatives u t (-,t) and f X xxxx( - , £) are from L°°([0, T]). By Gronwall's inequality, 
(14. Tl) yields w(t) = 0. Obviously, these estimates can be translated to the continuous 
dependence result in L 2 and hence in L\ oc . 

Other a priori bounds on solutions can be also derived in the lines of computations in 
[71 §§ 2,3] that lead to rather technical manipulations. The principal fact is the same as 
seen from (14.71) : differentiating a times in x equation (14. ip and setting v = D"u yields 
the equations with the same as in (14.51) principal part: 

(4.8) v t = -uv xxxxx + ... . 

Multiplying this by with ( being a cut-off function, and using various interpolation 
inequalities makes it possible to derive necessary a priori bounds and hence to observe 
the corresponding smoothing phenomenon for exponentially decaying initial data. 
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4.3. On local semigroup of smooth solutions of uniform NDEs and linear oper- 
ator theory. We recall that local C^-smoothing phenomena are known for third-order 
linear and fully nonlinear dispersive PDEs; see [H [TJ [22j [27] and earlier references therein. 
We claim that, having obtained a priori bounds, a smooth local solution can be con- 
structed by the iteration techniques as in [7J § 3] by using a standard scheme of iteration 
of the equivalent integral equation for spatial derivatives. We present further comments 
concerning other approaches to local existence, where we return to integral equations. 

We then need a detailed spectral theory of fifth-order operators such as 

(4.9) P 5 = a(x)D 5 x + b(x)D* + ..., xe(-L,L) (a(x)>±>0), 

with bounded coefficients. This theory is developed in Naimark's book [291 Ch. 2]. For 
regular boundary conditions (e.g., for periodic ones that are regular for any order and that 
suit us well), operators (14.91) admit a discrete spectrum {A^}, where the eigenvalues A& 
are all simple for all large k. 

It is crucial for further use of eigenfunction expansion techniques that the complete in 
L 2 subset of eigenfunctions {ipk} creates a Riesz basis, i.e., for any / G L 2 , 

(4.10) EI(/,^)| 2 <oo, where (/,V*> = 

and, for any {c^} G I 2 (i.e., ^ l c fc| 2 < oo), there exists a function / G L 2 such that 

(4.11) (fM = c k . 

Then there exists the unique set of "adjoint" generalized eigenfunctions {ipD (attributed 
to the"adjoint" operator Pg) being also a Riesz basis that is bi-orthonormal to {ipk}'- 

(4.12) {ip k ^i) = ki- 

Hence, for any / G L 2 , in the sense of the mean convergence, 

(4.13) f = Y J c k ijj k , with c k = (f,i>* k ). 

See further details in [291 § 5]. 

The eigenvalues of (14. 9 j) have the asymptotics 

(4.14) A fc ~ (±2vrA;i) 5 for all k » 1. 

In particular, it is known that P5 has compact resolvent that makes it possible to use it 
in the integral representation of the NDEs; cf. |7, § 3], where integral equations are used 
to construct a unique smooth solution of third-order NDEs. 

On the other hand, this means that P5— al for any a ^> 1 is not a sectorial operator that 
makes suspicious using advanced theory of analytic semigroups [SI ESI ES] , which is natural 
for even-order parabolic flows; see further discussion below. Analytic smoothing effect for 
higher-order dispersive equations were studied in [IS]. Concerning unique continuation 
and continuous dependence properties for dispersive equations, see [TTJ and references 
therein, and also [H] for various estimates. 
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For the linear dispersion equation with constant coefficients f !2. 13j) . the Cauchy problem 
with integrable data u (x) admits the unique solution 

(4.15) u(x,t) = b(x — -,t) * uo(-), 

where b(x,t) is the fundamental solution (12.141) . Analyticity of solutions in t (and x) can 
be associated with the rescaled operator 

(4.16) B 5 = -D 5 z + ±zD z + ll in L 2 p (R), p(z) = e°l*l 8/4 , 

where a > is a sufficiently small constant. Here, B5 in (I4.16P is the operator in (12.151) 
that generates the rescaled kernel F of the fundamental solution in (12.141) . Namely, using 
the same rescaling as in (12.141) . we set 

(4.17) u(x, t) — t~$v{y, t), y — x/ts, r = lnt, 
to get the rescaled PDE with the operator (14.161) . 

(4.18) v T = B 5 v, 

so that ( 14.15p . on Taylor's expansion of the kernel, yields 

(4.19) v(y, r) = j F(y - ze-/ 5 ) u (z) dz = Z^ F^(y)e-^^ j z k u (z) dz, 

(k) 

where the series converges uniformly on compact subsets thus defining an analytic solution. 
Then 1 14. 16ft has the real spectrum and the eigenf unctions (see [TSl § 9]) 

or(B 5 ) = {-i, Z = 0,1,2,...} and Mv) = ^ F {k Kv)> k >°- 

See a "parabolic" version of such a spectral theory in [12]. This suggests that B 5 — al 
is sectorial for a > (A = is simple), and this justifies the fact that (14.151) is an 
analytic (in t) flow. Let us mention again that analytic smoothing effects are known for 
higher-order dispersive equations with operators of principal type, [43]. 

Actually, this suggests to treat (14. ip . (14. 2 p by classic approach as in Da Prato-Grisvard 
[8] by linearizing about a sufficiently smooth m = u(to), to > 0, by setting u(t) = uo+v(t) 
giving the linearized equation 

(4.20) v t = A'(u )v + A(u )+g(v), t > t ; v(t o ) = 0, 

where g(v) is a quadratic perturbation. Using the good semigroup e A '( U0 ^, this makes it 
possible to study local regularity properties of the corresponding integral equation 

t 

(4.21) v(t) = J e A '( u °)(*- s ) (A(u ) + g(v(s))) ds. 

to 

Note that this smoothing approach demands a fast exponential decay of solutions v(x,t) 
as x — > 00, since one needs that v (-, t) G L 2 p ; cf. [2Z], where C^-smoothing for third-order 
NDEs was also established under the exponential decay. Equation (14.211) can be used to 
guarantee local existence of smooth solutions of a wide class of odd-order NDEs. 
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Thus, we state the following conclusion to be used later on: 



any sufficiently smooth solution u(x, t) of ( 14. ip . (14. 2ft at t — t , 
can be uniquely extended to some interval t £ (t ,t + v), v > 0. 

4.4. Smooth deformations and 5-entropy test for solutions with shocks. The 

situation dramatically changes if we want to treat solutions with shocks. Namely, it is 
known that even for the NDE-3 (jl.28p . the similarity formation mechanism of shocks 
shows nonuniqueness extension of solutions after a typical "gradient" catastrophe, |17j . 
Therefore, we dot a chance to get in an easy (or any) manner a uniqueness/entropy 
result for more complicated NDEs such as (II. 5p by using the 5-deformation (evolutionary 
smoothing) approach. However, we will continue using these fruitful ideas in order to 
develop a much weaker "5 -entropy test" for distinguishing various shock and rarefaction 
waves. 

Thus, given a small 5 > and a sufficiently small bounded continuous (and, possibly, 
compactly supported) solution u(x,t) of the Cauchy problem (14. ip . satisfying (14.21) . we 
construct its smooth 5 -deformation, aiming get smoothing in a small neighbourhood of 
bounded shocks, as follows. Note that we deal here with simple shock configurations 
(mainly, with 1-shock structures), and do not aim to cover more general shock geometry, 
which can be very complicated; especially since we do not know all types of simple single- 
point moving shocks. 

(i) We perform a smooth ^-deformation of initial data Uq(x) by introducing a suitable 
C 1 function u s(x) such that 

(4.23) J\u -u os \<5. 

If Mo is already sufficiently smooth, this step is abandoned (now and later on). By uu(x, t), 
we denote the unique local smooth solution of the Cauchy problem with data uqs, so that, 
by (14.221) . continuous function ui$(x,t) is defined on the maximal interval t G [t ,ti(S)), 
where we denote to — and ti(S) = Ai$. At this step, we are able to eliminate non- 
evolution (evolutionary unstable) shocks, which then create corresponding smooth rar- 
efaction waves. 

(ii) Since at t = A^, a shock-type discontinuity (or possibly infinitely many shocks) 
is supposed to occur, since otherwise we extend the continuous solution by (I4.22p . we 
perform another suitable <5- deformation of the "data" uis(x, Ai$) to get a unique contin- 
uous solution v,2s(x,t) on the maximal interval t 6 [ti(5), t 2 (^)), with t 2 {5) = Ais + A 2 s, 
etc. Here and in what follows, we always mean a "5-smoothing" performed in a small 
neighbourhood of occurring singularities only (shock waves). 

(k) With suitable choices of each 6- deformations of "data" at the moments t = tj(S), 
when Uj$(x,t) has a shock, there exists a £&(£) > 1 for some finite k = k(5), where 
k(5) — > +oo as 5 — > 0. It is easy to see that, for bounded solutions, k(5) is always finite. 
A contradictions is obtained while assuming that tj (5) — > i < 1 as j ; — > oo for arbitrarily 
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small 5 > meaning a kind of "complete blow-up" that was excluded by assumption of 
smallness of the data. 

This gives a global 5 -deformation in R x [0,1] of the solution u(x,t), which is the 
discontinuous orbit denoted by 

(4.24) u s (x, t) = {u j5 {x, t) for t G tj($)) t J = 1,2,..., k(6)}. 

One can see, this 5-deformation construction aims checking a kind of evolution stability 
of possible shock wave singularities and therefore, to exclude those that are not entropy 
and evolutionary generate smooth rarefaction waves. 

Finally, by an arbitrary smooth 5 -deformation, we will mean the function (14.241) con- 
structed by any sufficiently refined finite partition {tj(5)} of [0, 1], without reaching a 
shock of SL-type at some or all intermediate points t = tj(5). 

We next say that, given a solution u(x,t), it is stable relative smooth deformations, or 
simply 5-stable (5 eformation- stable) , if for any e > 0, there exists 5 = 5(e) > such that, 
for any finite 5-deformation of u given by (14.241) . 

(4.25) JJ\u-u 5 \<£. 

Recall that (I4.24p is an 5-orbit, and, in general, is not and cannot be aimed to represent 
a fixed solution in the limit 5 — > 0; see below. 

4.5. On 5-entropy solutions. Having checked that the local smooth solvability problem 
above is well-posed, we now present the corresponding definition that will be applied to 
particular weak solutions. Recall that the metric of convergence, L\ oc in present, for 
(11.281) was justified by a similarity analysis presented in Proposition 12.11 For other types 
of shocks and/or NDEs, the topology may be different. 

Thus, under the given hypotheses, a function u(x,t) is called a 5-entropy solution 
of the Cauchy problem (14.11) . if there exists a sequence of its smooth ^-deformations 
{u &k , k = 1, 2, ...}, where 5^ — * 0, which converges in L\ oc to u as k — > oo. 

This is slightly weaker (but equivalent) to the condition of 5-stability. 

4.6. ^-entropy test and on nonexistent uniqueness. Since, by obvious reasons, the 
5-deformation construction gets rid of non-evolution shocks (leading to non-singular rar- 
efaction waves), a first consequence of the construction is that it defines the 5- entropy 
test for solutions, which allows one, at least, to distinguish the true shocks from artificial 
rarefaction wave. 

In Section O we show that it is completely non-realistic to expect something essentially 
more from this construction in the direction of uniqueness and/or entropy-like selection 
of proper solutions. Though these expectations well-correspond to previous classic PDE 
entropy-like theories, these are too much excessive for higher-order models, where such a 
universal property is not achievable at all anymore. Even proving convergence for a fixed 
special 5-deformation is not easy at all. Thus, for particular cases, we will use the above 
notions with convergence along a subsequence of £'s to classify and distinguish shocks and 
rarefaction waves: 
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4.7. First easy conclusions of 5-entropy test. As a first application, we have: 



Proposition 4.1. Shocks of the type S-(x) are 5 -entropy for (4-1)- 

The result follows from the properties of similarity solutions ( 12. ip . which, by varying 
the blow-up time T i— > T + 5, can be used as their local smooth (^-deformations at any 
point t G [0, 1). 

Proposition 4.2. Shocks of the type S+(x) are not 5 -entropy for ( 4-1 )- 

Indeed, taking initial data S + (x) and constructing its smooth ^-deformation via the 
self-similar solution (13.21) with shifting t i— > t + 5, we obtain the global 5-deformation 
{u s = u+(x, t + 5)}, which goes away from S + . 

Thus, the idea of smooth 5-deformations allows us to distinguish basic 5-entropy and 
non-entropy shocks without any use of mathematical manipulations associated with stan- 
dard entropy inequalities, which, indeed, are illusive for higher-order NDEs; cf. [17j. 

5. On nonuniqueness after shock formation 

Here we mainly follow the ideas from [T7j applied to the NDE-3 (11.281) . so we will omit 
some technical data and present more convincing analytic and numerical results (numerics 
are essential for occurred hard 5D dynamical systems). Without loss of generality, we 
always deal with the divergence NDE-5 (|1.5p . 

5.1. Main strategy towards nonunique continuation. We begin with the study of 
new shock patterns, which are induced by other similarity solutions of (II. 5p : 

(5.1) u.(x,t) = (-t) a f(y), y = j^, P = ^, where a G (0, |) and 



(5.2) 



-(ff') W -Pf'y + af = in R_, /(0) = f'(0) = / (4) (0)=0, 
f(v)=C \y\?(l + o(l)) as y^-oo, C > 0. 



In this section, in order to match the key results in [T7j, in (12. ip and later on, we change 
the variables {g,z} i— > {/, y}. In the next Section [6j we return to the original notation. 

The anti-symmetry conditions in ( 15. 2p allow to extend the solution to the positive 
semi-axis {y > 0} by —f(—y) to get a global pattern. 

Obviously, the solutions (I2.ip . which are suitable for Riemann's problems, correspond 
to the simple case a = in (15.11) . It is easy to see that, for positive a, the asymptotics in 
( 15. 2p ensures getting first gradient blow-up at x = as t —>■ CT, as a weak discontinuity, 
where the final time profile remains locally bounded and continuous: 

(5.3) u_(z,0~) = 

where Co > is an arbitrary constant. Note that the standard "gradient catastrophe r 
u x (0, 0~) = — oo, then occurs in the range, which we will deal within, 

(5.4) | < 1 provided that a < |. 
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Thus, the wave braking (or "overturning") begins at t — 0, and next we show that it is 
performed again in a self-similar manner and is described by similarity solutions 

(5.5) u + (x,t)=t a F(y), y = $, (3 = i±2, where 

(-(FF>)^+f3F'y-aF = in R_, 

{ F(0) = F >0, F(y) = C |y|t(l + o(l)) as y -> -oo, 

where the constant Co > is fixed by blow-up data (15.31) . The asymptotic behaviour 
as y — > — oo in (15. 6p guarantees the continuity of the global discontinuous pattern (with 
F(—y) = —F(y)) at the singularity blow-up instant t = 0, so that 

(5.7) u_(ar,0 - ) = u + {x,0 + ) in R. 

Then any suitable couple {/, F} defines a global solution u±(x, t), which is continuous 
at t = 0, and then it is called an extension pair. It was shown in [TTj that, for typical 
NDEs-3, the pair is not uniquely determined and there exist infinitely many shock-type 
extensions of the solution after blow-up at t — 0. We are going to describe a similar 
non-uniqueness phenomenon for the NDEs-5 such as (II. 5p . 

A first immediate consequence of our similarity blow- up/extension analysis is as follows: 

(5.8) in the CP, formation of shocks for the NDE (11.51) can lead to nonuniqueness. 

The second conclusion is more subtle and is based on the fact that, for some initial data 
at t — 0, the solution set for t > does not contain any "minimal", "maximal", or 
"extremal" points in any reasonable sense, which might play a role of a unique "entropy" 
one chosen by introducing a hypothetical entropy inequalities, conditions, or otherwise. If 
this is true for the whole set of such weak solutions of (11.51) with initial data (15.31) . then, 
for the Cauchy problem, 

(5.9) there exists no general "entropy mechanism" to choose a unique solution. 

Actually, overall, (I5.8P and (15.91) show that the problem of uniqueness of weak solutions 
for the NDEs such as (11.51) cannot be solved in principal. On the other hand, in a FBP 
setting by adding an extra suitable condition on shock lines, the problem might be well- 
posed with a unique solution, though proofs can be very difficult. We refer again to more 
detailed discussion of these issues for the NDE-3 (11.281) in [17]. Though we must admit, 
that for the NDE-5 (II. 5p . which induces 5D dynamical systems for the similarity profiles 
(and hence 5D phase spaces), those nonuniqueness and non-entropy conclusions are not 
that clear and some of their aspects do unavoidably remain questionable and open. 

Hence, the non-uniqueness in the CP is a non-removable issue of PDE theory for higher- 
order degenerate nonlinear odd-order equations (and possibly not only those). The non- 
uniqueness of solutions of (11.51) has a pure dimensional nature associated with the struc- 
ture of the 5D phase space of the ODE (15. 6p . 
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5.2. Gradient blow-up similarity solutions. Consider the blow-up ODE problem 
(15. 2p . which is a difficult one, with a 5D phase space. Note that, by invariant scaling 
(12.211) . it can be reduced to a 4th-order ODE with a too complicated nonlinear operator 
composed from too many polynomial terms, so we do not rely on that and work in the 
original phase space. Therefore, some more delicate issues on, say, uniqueness of certain 
orbits, become very difficult or even remain open, though some more robust properties 
can be detected rigorously. We will also use numerical methods for illustrating and even 
justifying some of our conclusions. For the fifth-order equations such as (I5.2p . this and 
further numerical constructions are performed by the MatLab with the standard ode45 
solver therein. 

Let us describe the necessary properties of orbits {f(y)} we are interested in. Firstly, 
it follows from the conditions in (15.21) that, for y « 0~, 

(5.10) the set of proper orbits is 2D parameterized by fi = f'(0) < and / 3 = f"'(0). 
Secondly and on the other hand, the necessary behaviour at infinity is as follows: 

(5.11) /(2/) = CoM^(l + o(l)) as y^-oo = f), 

where Co > is an arbitrary constant by scaling (12.211) . It is key to derive the whole 4D 
bundle of solutions satisfying (15. lip . This is done by the linearization as y — > — oo: 

f(y) = fo(y)+Y(y), where f (y) = C (-y)f 

=> -C ((-y)^Y)^ + (3Y'(-y) + aY + l(f^y))^ + ... = 0. 

By WKBJ-type asymptotic techniques in ODE theory, solutions of (I5.12p have a standard 
exponential form with the characteristic equation: 

(5.13) Y{y) ~ e a ^)\ 7 = 1 + ± (l - f ) > 1 =► C ( 7 a) 4 = /?, 
which has three roots with Rea^ < 0. Hence, we conclude that: as y — > — oo, 

(5.14) the bundle (15.111) is four-dimensional. 

The behaviour (15.151) with the bundle (I5.14p gives the desired asymptotics: by (15. lip , we 
have the gradient blow-up behaviour at a single point: for any fixed x < 0, as t — ► CT, 
where y = x/(—tf — > — oo, uniformly on compact subsets, 

(5.15) u.{x,t) = (-t) a f(y) = (-t) a C \-^\%(l + o(l)) - C \x\^, 

Let us explain some other crucial properties of the phase space, now meaning "bad 
bundles" of orbits. First, these are the fast growing solutions according to the explicit 
solution 

(5.16) Uy) = -jg^ > for y « -1. 
Analogously to (15 . 1 2j) . we compute the whole bundle about ( 15.161) : 

(5.17) f(y) = f*(y)+Y(y) =}► TE ^ -(y 5 Y)^-f3Y'y + aY + ... = 0. 
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Blow-up f(y): function h (m) for a=1/9 with five real negative roots 



Blow-up f(y): function h (m) for a=1/9, 1/19, and 3/17 (five real roots) 




(a) a = \ 




(b) various a 



FIGURE 9. The polynomial h a (m) in (I5.18|) for various a £ (0, |): five negative roots. 



This Euler's equation has the following solutions with the characteristic polynomial: 

(5.18) Y(y) = y m =^ h a (m) = ("+i)(^+2)(m+3)(m+4)(m+ 5 ) _ (3m + a = Q _ 

One root m = —5 is obvious that gives the solution (15.161) . It turns our that this algebraic 
equation has precisely five negative real roots for a from the range (15.41) . as Figure [9] shows. 
Actually, (b) explains that the graphs are rather slightly dependent on a. Thus: 

(5.19) the bundle about (15.161) is five-dimensional. 

Second, there exists a bundle of positive solutions vanishing at some finite y — > y$ < 
with the behaviour (this bundle occurs from both sides, as y — > y$ to be also used) 

(5.20) /i(y)= VbFjfol (l + o(l)), A > 0, 

is 4D, which also can be shown by linearization about (15.201) . Indeed, the linearized 
operator contains the leading term 

(5.21) -A 2 (^\y~yV\Y)^ + ... = Y(y) ~ \y - y \l, \y - y \K \y - y \h 
which together with the parameter y < yields 

(5.22) the bundle about (15.201) is four- dimensional. 

Thus, ( 15.101) . (I5.14p . (I5.19p . and (15.221) prescribe key aspects of the 5D phase space we 
are dealing with. To get a global orbit {f(y),y G K_} as a connection of the proper 
bundles (15.10p and ( 15. 14ft . it is natural to follow the strategy of "shooting from below" by 
avoiding the bundle ( 15.20D . ( 15.221) . i.e., using the parameters /i^ in ( 15.101) . to obtain 

(5.23) y = -oo. 

It is not difficult to see that this profile f(y) will belong to the bundle (15.141) . The proof 
of such a 2D shooting strategy can be done by standard arguments. By scaling (I2.2ip . we 
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Typical shooting of blow-up profile f(y) fora=1/9: f Q = f 2 =f 4 =0, ^=-1 




FIGURE 10. The shooting strategy of a blow-up similarity profile f(y) for a = ^, 
with data /(0) = /"(0) = = and /'(0) = 1; the shooting parameter is 
/'"(0) = 0.0718040128557... . 

always can reduce the problem to a ID shooting: 

(5.24) / 1 = /'(0) = -l and / 3 = /'"(0) is a parameter. 

By the above asymptotic analysis of the 5D phase space, it follows that: 

(i) for /3 — 1 the orbit belongs to the bundle about (15.161) . and 

(ii) for / 3 3> 1, the orbit vanishes in finite y along (15.201) . 

Hence, by continuous dependence, we obtain a solution f(y) by the min-max principle 
(plus some usual technical details that can be omitted). Before stating the result we can 
prove, in Figure [TOl we present how we are going to justify existence of a proper blow-up 
shock profile f(y). Thus, we fix the above speculations as: 



Proposition 5.1. (i) In the range ( 5.J$ , the problem ( I5.#p admits a shock wave profile 



f(y); and (ii) this f(y) is unique up to scaling {\2.21\i and is positive for y < 0. 

Two results (i) and (ii) are combined for convenience, and now we declare that (ii) is 
an open problem. In [17], for the NDE-3 (11.281) . the phase space is 3D and the proof is 
available. 

In fact, this is a rather typical result for higher-order dynamical systems. E.g., we refer 
to a similar and not less complicated study of a 4th-order ODE [21] , where existence and 
uniqueness of a positive solution of the radial bi-harmonic equation with source: 

(5.25) A 2 r u = u p for r = |x| > 0, w(0) = 1, w'(0) = u'"(0) = 0, u(oo) = 0, 

was proved in the supercritical Sobolev range p > psob = j^zz? N > A. Here, analogously, 
there exists a single shooting parameter being the second derivative at the origin U2 = 
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FIGURE 11. Shooting the blow-up profile f(y) for a = |; here /(0) = 10. 

u"(0); the value tt = u(0) = 1 is fixed by a scaling symmetry. Proving uniqueness of such 
a solution in [21] is not easy and lead to essential technicalities, which the attentive Reader 
can consult in case of necessity. Fortunately, we are not interested in any uniqueness of 
such kind. Instead of the global behaviour such as ( 15.161) . the equation ( 15.251) admits the 
blow-up one governed by the principal operator u^' + ... = u p (u — > +oo). The solutions 
vanishing at finite point otherwise can be treated as in the family (I). 

We next use more advanced and enhanced numerical methods towards existence (and 
uniqueness-positivity, see (ii)) of f(y). Figure [TT1 shows the shooting from y = 0~ for 

(5.26) a = \ | = i 

This again illustrates the actual strategy in proving Proposition 15.11 Note that here, as 
an illustration of another important evolution phenomenon, we solve the problem with 

(5.27) /(0) = / = 10 =^ MO,t)] = 2/o(-t) a ^0 as * - 0" 

so that this similarity solution describes collapse of a shock wave. 

Next, Figure [121 shows truly blow-up profiles, with /(0) = 0, constructed by a different 
method (via the solver bvp4c) for convenient values a = |, ^, and Note the clear 
oscillatory behaviour of such patterns that is induces by complex roots of the characteristic 
equation (15.131) . 

Stationary solutions with a "weak shock" . The ODE in (pT2l) and hence the 
PDE (ll.5p admit a number of simple continuous stationary solutions. E.g., consider 

(5.28) a = §, | = | : f{y) = y^\y\ sign y and u(x, t) = a/|x[ signx. 
We will show that such weak shocks also lead to nonuniqueness. 
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Shock formation blow-up similarity profiles f(y) fora>0: f Q =f 2 =f 4 =0 

25 I 1 1 1 1 1 1 




y 



FIGURE 12. The odd blow-up similarity profiles f(y) in M_ with a = g, yg, and ^. 



5.3. On nonuniqueness of similarity extensions beyond blow-up. A discontinuous 
shock wave extension of blow-up solutions (15.11) . (I5.2p is done by using the global ones 
(15. 5p . (15. 6p . Actually, this leads to watching a whole 5D family of solutions parameterized 
by their Cauchy values at the origin: 

(5.29) F(0) = F >0, F'(0) = F 1 <0, F"(0) = F 2 , F"'(0) = F 3 , F (4) (0) = F A . 

Thus, unlike (I5.10p . the proper bundle in (I5.29j) is 5D. Note that at minus infinity, the 
solution must have the form 

(5.30) F{y) = C \y\^{l + o{l)) as y -> -oo (C > 0). 

As above, the 5D phase space for the ODE in (I5.6P has two stable "bad" bundles: 

(I) Positive solutions with "singular extinction" in finite y, where F(y) —>■ as y —>■ y$ < 
0. This is an unavoidable singularity following from the degeneracy of the equations with 
the principal term FF^ leading to the singular potential ~ As in (15.211) . this bundle 
is 4D, and 

(II) Negative solutions with the fast growth (cf. (I5.16P ): 

(5.31) F*(y) = (1 + o(l)) ^ -oo as y ->• -oo. 

The characteristic polynomial is the same as in (I5.18p . so that the bundle is 5D; cf. (I5.19p . 

Both sets of such solutions are open by the standard continuous dependence of solutions 
of ODEs on parameters. The whole bundle of solutions satisfying (15.111) is obtained by 
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Nonexistence of global profile F(y) for a=1/9: F Q =1 , F =1, F =F =0 




Nonexistence of global profile F(y) for a=1/9: F Q =1, F 1= F 2 =F 3 =F 4 




(a) Fx = -1, F 2 = F 3 = 



(b) F 1 ^F 2 ^F 3 = F 4 



FIGURE 13. Examples of nonexistence of F(y) via shooting from y = 



linearization as y 
(5.32) 



0. 



-co in ( 15. 6ft : 

f(y) = F (y)+Y(y), where F (y) = G)(-# 

-c ((-y)* *1 (5) - /n"(-y) - or + J (F 2 ( y ))( 5 ) + ... 

The WKBJ method now leads to a different characteristic equation: 

(5.33) ~ e«~W, 7 = 1 + J (l - |) > 1 =s> C ( 7 a) 4 = 

so that there exist just two complex conjugate roots with Re< 0, and hence, unlike (I5.14p . 

(5.34) the bundle (15.111) of global orbits {F(y)} is three-dimensional. 

However, the geometry of the whole phase space changes dramatically in comparison 
with the blow-up cases, so that the standard shooting of positive global profiles F(y) 
yields no encouraging results. We refer to Figure [T3J which illustrate typical results of a 
standard shooting. In Figure [T3J we show two other examples showing nonexistence of 
such an F(y). In (a), this is done by shooting from from the left from the bundle ( I5.30p . 
and in (bO we have used the bvp4c solver to get convergence again to profiles with a 
singularity as in (15.201) . 

Thus, surprisingly, we arrive at a clear nonexistence of the extension similarity problem 
(I5.6p . Note that for the NDE-3 (11.281) . such a nonunique profile F(y) does exist [17], thus 
ensuring the nonuniqueness (and non-entropy) in the similarity-ODE framework for the 
problem. We pose two related questions: 

Open Problems [5ll (i) To prove that the global similarity extension problem ( \5.6§ does 
not have a solution F(y) for any Co > 0, and hence 

(ii) To describe non- self- similar formation of (nonunique) shock waves from data (\5.3§ . 
Nonuniqueness. The justified nonuniqueness is achieved for the values of parameters 

no) 



F(0) = F < 



and 
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Fi < 0, 



Nonexistence of global F(y) (ora=1/9: shooting from the left 



Nonexistence of global F(y) fora=1/9: F =10, F =-0.1, 0, +0.1, F =0 





(a) shooting from the left by ode45 



(b) by bvp4c 



FIGURE 14. Examples of nonexistence of F(y) for a = |: by shooting from the 
left (a) and by iteration techniques (b). 

as shown in Figure [Tj3 The proof of existence of such profiles F is based on the same 
geometric arguments as that of Proposition 12.21 (with the evident change of the geometry 
of the phase space). These two different profiles show a nonunique way to get solutions 
with the initial data (Co = 1 by scaling) 



u (x) 



\x\0signx in 



This is another type of nonuniqueness in the Cauchy problem for (11.51) . showing the 
nonunique way of formation of shocks from weak discontinuities, including the stationary 
ones as in (15. 28ft . 

5.4. More on nonuniqueness and well-posedness of FBPs. The non-uniqueness 
(15. 8p in the Cauchy problem (11.51) . (15.31) is: any F(y) yields the self-similar continuation 
(15. 5p . with the behaviour of the jump at x = (profiles F(y) as in Figure [T5|) 

(5.35) 

In the similarity ODE representation, this nonuniqueness has a pure dimensional origin 
associated with the dimension of the good and bad asymptotic bundles of the 5D phase 
spaces of both blow-up and global equations. Since these shocks are stationary, the cor- 
responding Rankine-Hugoniot (R-H) condition on the speed A of the shock propagation: 



[u + (x,t)]\ x=0 = -(u + (0 + ,t)-u + (0-,t)) = 2F t a <0 for t > 0. 



(5.36) 



A 



[(mii)i 



\x=Q 







is valid by anti-symmetry. As usual, (I5.36P is obtained by integration of the equation 
(11.281) in a small neighbourhood of the shock. The R-H condition does not assume any 
novelty and is a corollary of integrating the PDE about the line of discontinuity. 

Moreover, the R-H condition (I5.36P also indicates another origin of nonuniqueness: a 
symmetry breaking. Indeed, the solution for t > is not obliged to be an odd function of 
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Shooting global profile F(y) for o=1/9: F =-1, F=F =-0.115526..., F =F =0 



Global profile F(y) fora=1/9: F =-1, F =0, F =-0.16648, F =F =0 




(a) F'(0) =F"'(0) 



-0.115526. 




(b) F'(0) = 0, F"(0) 



FIGURE 15. Shooting a proper solution F(y) of (15. 6|) for a 



-0.16648. 



with data 



F(0) = -1, F A = 0, and F'(0) = F"'(0) 
F'(0) = 0, F 2 = -0.16648... , F 3 = 0. 



-0.115526..., F 2 = (a), and 



x, so the self similar solution ( 15.5ft for x < and x > can be defined using ten different 
parameters {i^o , F^}, and the only extra condition one needs is the R H one: 

(5.37) [(F F')"'\ = 0, i.e., F ~F^ + AF{ F^ + 3(F 2 ~) 2 = + 4F+F 3 + + 3(F+) 2 . 

This algebraic equations with ten unknowns admit many other solutions rather than the 
obvious anti-symmetric one: 

F - = -F+, FT = F+, F 2 ~ = -Ft, F 3 - = F 3 + , and F 4 " = 

Finally, we note that the uniqueness can be restored by posing specially designed con- 
ditions on moving shocks, which, overall guarantee the unique solvability of the algebraic 
equation in (I5.37P and hence the unique continuation of the solution beyond blow-up. 
This construction is analytically similar to that for the NDEs-3 (11.281) in [T7] . 

6. Shocks for an NDE obeying the Cauchy-Kovalevskaya theorem 

In this short section, we touch the problem of formation of shocks for NDEs that are 
higher-order in time. Instead of studying the PDEs such as (cf. [16l [19] ) 



(6.1) 



u u = ~{uu x ) 3 



Uttt = -{uUx)z 



etc. 



we consider the fifth-order in time NDE (ll.lOp . which exhibits certain simple and, at the 
same time, exceptional properties. Writing it for W = (u,v,w,g, h) T as 



(6.2) 



u t -- 
v t ~- 
w t 
9t -- 
h t 



v x , 
w x , 

- 9x, 
h 

■ UUn 



or Wt = AW T , with A 
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(j 1.1 Op becomes a first-order system with the characteristic equation for eigenvalues 

A 5 -u = 0, 

so it is not hyperbolic for 

6.1. Evolution formation of shocks. For (11.101) . the blow-up similarity solution is 

(6.3) U-(x, t) — g(z), z — x/(—t), where 

(gg')W = (z 5 g')^ = UOg'z + 2AQg"z 2 + 120g"'z 3 

(6.4) 

+20g i - 4 h 4 + g^z 5 in R, /(^oo) = ±1. 
Integrating (16 .4p four times yields 

(6.5) gg' = z 5 g' + Az + Bz 3 , with constants A = (g'(0)) 2 > 0, B = § g'(0)g'"(0), 
so that the necessary similarity profile g(z) solves the first-order ODE 

(e-e) s = ^sr ■ 

By the phase-plane analysis of (16.61) with A > and B = 0, we easily get the following: 

Proposition 6.1. The problem {6Jj_) admits a solution g(z) satisfying the anti-symmetry 
conditions ( \2. 9\\ that is positive for z < 0, monotone decreasing, and is real analytic. 

Actually, involving the second parameter B > yields that there exist infinitely many 
shock similarity profiles. The boldface profile g(z) in Figure [TBI (by (16.31) . it gives S-(x) 
as t — > 0~) is non-oscillatory about ±1, with the following algebraic rate of convergence 
to the equilibrium as z —>■ — oo: 

T + 3^ + ... for B = 0, 
1 + £ + ... for B >0. 

Z 

Note that the fundamental solutions of the corresponding linear PDE 
(6.7) 

is not oscillatory — > dboo. This has the form 

b(x,t) =t 3 F(y), y = x/t, so that b(x,0) = ... = b m (x,0) = 0, fo ttt t(rc, 0) = 8(x). 

The linear equation ( 16. 7ft exhibits some features of finite propagation via TWs, since 

u(x,t) = f(x - Xt) =}► -A 5 / {5) = / {5) , i-e., A = -1 

and the profile f(y) disappears from. This is similar to some canonical equations of 
mathematical physics such as 

Ut = u x (dispersion, A = — 1) and u u = u xx (wave equation, A = ±1). 

The blow-up solution (16.31) gives in the limit t — > 0~ the shock S- (x) , and (11.71) holds. 

Since (11.101) has the same symmetry (13. lj as (11.281) . similarity solutions (16.31) . with 
— t t— > t and g(z) h- > g(—z) according to (13. 5p . also give the rarefaction waves for S+(x) as 
well as other types of collapse of initial non-entropy discontinuities. 
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FIGURE 16. The shock similarity profile satisfying (|6.4p , 

6.2. Analytic 5-deformations by Cauchy— Kovalevskaya theorem. The great ad- 
vantage of the equation (ll.lOp is that it is in the normal form, so it obeys the Cauchy- 
Kovalevskaya theorem [151 p. 387]. Hence, for any analytic initial data u(x,0), u t (x,0), 
u tt (x,0), u m (x,0), and u U tt(%, 0), there exists a unique local in time analytic solution 
u(x, t). Thus, (jl.lOp generates a local semigroup of analytic solutions, and this makes it 
easier to deal with smooth 5-deformations that are chosen to be analytic. This defines 
a special analytic 5-entropy test for shock/rarefaction waves. On the other hand, such 
nonlinear PDEs can admit other (say, weak) solutions that are not analytic. Actually, 
Proposition 16.11 shows that the shock S-(x) is a 5-entropy solution of (11.101) . which is 
obtained by finite-time blow-up as t —>■ 0~ from the analytic similarity solution (16.31) . 

6.3. On formation of single-point shocks and extension nonuniqueness. Similar 
to the analysis in Section [5J, for the model (jl.lOp (and (16.11) ). these assume studying 
extension similarity pairs {/, F} induced by the easy derived analogies of the blow-up (15.21) 
and global (I5.6P 5D dynamical systems. These are very difficult, so that checking possible 
nonuniqueness and non-entropy of such flows with strong and weak shocks becomes a 
hard open problem (some auxiliary analytic steps towards nonuniqueness are doable). 

7. NONNEGATIVE COMPACTONS OF FIFTH-ORDER NDES ARE NOT ROBUST 

We begin with an easier explicit example of compactons for a third-order NDE. 

7.1. Third-order NDEs: 5-entropy compactons. Compactons as compactly sup- 
ported TW solutions of the A(2, 2) equation (11.221) were introduced in 1993, [39], as 

(7.1) u c (x,t) = f c (y), y = x + t => f c : f=(f)" + f. 
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Integrating yields the following explicit compacton profile: 



(7-2) f c (y) 



4 „„ n 2 



COS' 



(!) for \y\<2n, 



for \y\ > 2n. 

The corresponding compacton (17. ip . (17.21) is a 5-entropy solution, i.e., can be constructed 
by smooth (and moreover analytic) approximation via strictly positive solutions, |16j . 

It is curious that the same compactly supported blow-up patterns occur in the com- 
bustion problem for the related react ion- diffusion parabolic equation 

(7.3) u t = {u 2 ) xx + u 2 . 

Then the standing- wave blow-up (as t —* T~) solution of S-regime leads to the same ODE: 

(7.4) u s (x,t) = (T-t)- 1 f(x) => f = {f)" + f. 

This yields the Zmitrenko-Kurdyumov blow-up localized solution, which has been known 
since 1975; see more historical details in [201 § 4.2]. 

7.2. Examples of C 3 -smooth nonnegative compacton for higher-order NDEs. 

Such an example was given in [101 p. 4734]. Following [20, p. 189], we construct this 
explicit solution as follows. The operator ~F§(u) of the quinic NDE 

(7.5) u t = F 5 (u) = (u 2 ) xxxxx + 25(u 2 ) xxx + UA(u 2 ) x , 
is shown to preserve the 5D invariant subspace 

(7.6) W 5 = C{1, cosx, sinx, cos2x, sin2x}, 

i.e., F 5 (W5) C W5. Therefore, (17.51) restricted to the invariant subspace W$ is a 5D 
dynamical system for the expansion coefficients of the solution 

u(x,t) = Ci(t) + C 2 (t) cosx + C 3 (t) smx + C 4 (t) cos2x + C 5 (t) sin2x. 

Solving this yields the explicit compacton TW 

lib cos4 (I) for \y\ < 

for \y\ > 71. 

This C x solution can be attributed to the Cauchy problem for (17. 5ft . 

The above invariant subspace analysis applies also to the 7th-order PDE 

(7.8) u t = F 7 (u) = D 7 x (u 2 ) + (3Dl{u 2 ) + 7 (^W + v{u 2 ) x - 
Here F 7 admits W§, if 

= 25, 7 = 144, and u = 0. 
Moreover [201 P- 190], the only operator F 7 in (17.81) preserving the 7D subspace 

(7.9) W7 = C{1, cosx, sinx, cos2x, sin2x, cos3x, sin3x} 
is in the following NDE-7: 

(7.10) u t = F 7 (u) = D 7 x (u 2 ) + 77D 5 x (u 2 ) + lS76{u 2 ) xxx + 14400(u 2 ),. 

This makes it possible to reduce (I7.10P on W7 to a complicated dynamical system. 
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(7.7) u c (x,t) = f c (x + t), where f c (y) 



7.3. Why nonnegative compactons for fifth-order NDEs are not robust. As 

usual, by robustness of such solutions we mean that these are stable with respect to small 
perturbations of the parameters entering the NDE and the corresponding ODEs. In other 
words, NDEs admitting such non-negative "heteroclinic" orbits 0^0 are not structurally 
stable in a natural sense. This reminds the classic Andronov-Pontriagin-Peixoto theorem, 
where one of the four conditions for the structural stability of dynamical systems in M 2 
reads as follows [32[ p. 301]: 

(7.11) "(ii) there are no trajectories connecting saddle points... ." 

Actually, nonnegative compactons are special homoclinics of the origin, and we will show 
that the nature of their non-robustness is in the fact that these represent a stable-unstable 
manifold of the origin consisting of a single orbit. Therefore, in consistency with (17. lip , 
the origin is indeed a saddle in M 4 . 

In order to illustrate the lack of such a robustness, consider the NDE (17.51) . where, on 
integration, we obtain the following ODE: 

(7.12) u c (x,t) = f c (x + t) => f c : 2/ = (/ 2 ) (4) + -, 

where we omit the lower-order terms. Looking for the compacton profile f > 0, we set 
f 2 = F to get 



(7.13) F (4) = 2VF + ... for y > 0, F'(0) = F'"(0) = 0. 

As usual, we look for a symmetric F(y) by putting two symmetry conditions at the origin. 

Let y = yo > be the interface point of F(y). Then, looking for the expansion as 
y — > y^ in the form 

(7.14) F(y) = J ±?(y -y) a + e(y), with e(y) = o((y - y) 8 ), 
we obtain Euler's equation for the perturbation e(y), 

(7-15) i- (yo-y) 4 e^-e = 0. 

Hence, s(y) = (y — y) m , with the characteristic equation 

(7.16) m(m - l)(m - 2)(m - 3) - 840 = m x = -4, m 2 , 3 = 3±i ^ IIT , m 4 = 7. 

Hence, Rem ; < 8, and, in other words, (17.151) does not admit any nontrivial solution 
satisfying the condition in (17.141) : see further comments in [20l p. 142]. In fact, it is easy 
to see that ( 17.151) with e = is the unique positive smooth solution of F^ = 2yF. 

(7.17) Thus: the asymptotic bundle of solutions (17.141) is ID, 

where the only parameter is the position of the interface yo > 0. Obviously, as a typical 
property, this ID bundle is not sufficient to satisfy (by shooting) two conditions at the 
origin in (17.131) . so such TW profiles F(y) > are nonexistent for almost all NDEs like 
that. In other words, the condition of positivity of the solution, 



(7.18) to look for a nontrivial solution F > for the ODE in (17. 131) 
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creates a free-boundary "obstacle" problem that, in general, is inconsistent. Skipping 
the obstacle condition f 1 7 . 1 8 j) will return such ODEs (or elliptic equation), with a special 
extension, into the consistent variety, as we will illustrate below. 

Thus, nonnegative TW compactons are not generic (robust) solutions of (2m + l)th- 
order quadratic NDEs with m = 2, and also for larger m's, where (17.171) remains valid. 

7.4. Nonnegative compactons are robust for third-order NDEs only. The third- 
order case m — 1, i.e., NDEs such as (11.221) . is the only one where propagation of pertur- 
bations via nonnegative TW compactons is stable with respect to small perturbation of 
the parameters (and nonlinearities) of equations. Mathematically speaking, then the ID 
bundle in (17.171) perfectly matches with the single symmetry condition at the origin, 

F" = 2VF + ... and F'(0) = 0. 

7.5. Compactons of changing sign are robust. As a typical example, we consider 
the perturbed version (II. lip of the NDE-(1,4) (II. 5p . As we have mentioned, this is 
is written for solutions of changing sign, since nonnegative compactons do not exist in 
general. Looking for the TW compacton fl 7 . 1 2 [) yields the ODE 

(7.19) / = 4(|/|/) (4) + !l/|/ =► F^=F-2\F\^F (F=\f\f). 

Such ODEs with non-Lipschitz nonlinearities are known to admit countable sets of com- 
pactly supported solutions, which are studied by a combination of Lusternik-Schnirel'man 
and Pohozaev's fibering theory; see [18] . 

In Figure [HI we present the first TW compacton patterns (the boldface line) and the 
second one that is essentially non-monotone. These look like standard compacton profiles 
but careful analysis of the behaviour near the finite interface at y = yo shows that F(y) 
changes sign infinitely many times according to the asymptotics 

(7.20) F{y) = (y - y) 8 (p(s + s ) + s = ln{y - y). 

Here, the oscillatory component (f(s) is a periodic solution of a certain nonlinear ODE 
and s is an arbitrary phase shift; see [20, § 4.3] for further details. Thus, unlike (17. 17j) . 

(7.21) the asymptotic bundle of solutions (17.201) is 2D (parameters are yo and so), 

and this is enough to match two symmetry boundary conditions given in (17.131) . Such 
a robust solvability is confirmed by variational techniques that apply to rather arbitrary 
equations such as in (17.191) with similar singular non-Lipschitz nonlinearities. 

Regardless the existence of such sufficiently smooth compacton solutions, it is worth 
recalling again that, for the NDE (11.111) . as well as (I2.3ip and (14.31) both containing 
monotone nonlinearities, the generic behaviour, for other data, can include formation of 
shocks in finite time, with the local similarity mechanism as in Section [2j 
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Oscillatory compacton TW profiles for u =-(uu ) +uu , X=-1 
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FIGURE 17. First two compacton TW profiles F(y) satisfying the ODE in ([7719]) . 
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